

import gmpy2
import libnum
import sympy
n= 152362584260172832299347742092288781523553077595415336181107340981588793827198678351439908568365900532506619219107872595952010917469583778751598317914867718261964335080084208817370990836785169395479387595153363093720954122478858717743030228933613354156557566800833234418707803592840802973241337796008470981341
c= 19401395531981752916451282817866974743361193553919716050068261633260025385637451158918896346789157347763593932302429285603509319396908014135940103051580842310569537091180635065419629629004607071254512910448784377678421676801718187196005282018075668909717717835822216168107432901556029987260147722263316337224
p_q= 24761248958744201462825446279928534556231843440178848260305853629617417115589455163519829983799059130246977131703814834484201858080083363998431077557014558
e=65537
p = sympy.symbols('p')
q= sympy.symbols('q')
f1=p+q-p_q
f2=p*q-n
p,q=sympy.solve([f1,f2],[p,q],domain='ZZ')

print(p[0])
print(p[1])
phi=(p[0]-1)*(p[1]-1)
d=gmpy2.invert(e,int(phi))
m=pow(c,d,n)
print(libnum.n2s(int(m)))


